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Abstract 

Soft supersymmetry breaking terms involving heavy singlet sneutrinos provide new sources of 
lepton number violation and of CP violation. In addition to the CP violation in mixing, investigated 
previously, we find that 'soft leptogenesis' can be generated by CP violation in decay and in the 
interference of mixing and decay. These additional ways to leptogenesis can be significant for a 
singlet neutrino Majorana mass that is not much larger than the supersymmetry breaking scale, 
M < 10 2 msusY- In contrast to CP violation in mixing, for some of these new contributions the 
sneutrino oscillation rate can be much faster than the decay rate, so that the bilinear scalar term 
need not be smaller than its natural scale. 
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I. INTRODUCTION 



The evidence for neutrino masses at a scale ~ 10 2 eV makes a convincing case for 
the seesaw mechanism ODD : The existence of singlet neutrinos with Majorana masses 
and with Yukawa couplings to active neutrinos becomes very plausible. Thephysics of 
these heavy neutrinos can provide all the necessary ingredients for baryogenesis j^j: B — L is 
violated by the Majorana masses, CP is likely to be violated in the neutrino Yukawa couplings 
and, for small enough Yukawa couplings, the heavy neutrinos decay out of equilibrium. Thus, 
leptogenesis jj], the dynamical generation of lepton asymmetry through the decays of heavy 
singlet Majorana neutrinos, becomes an attractive solution to the puzzle of the baryon 
asymmetry. 

The seesaw mechanism introduces a new scale, M, the mass scale of the singlet neutrinos. 
Since this scale must be much higher than the electroweak breaking scale, M 3> Aew> a 
huge amount of fine-tuning is required within the framework of the Standard Model extended 
to include singlet neutrinos (SM+N) to keep the low Higgs mass. This situation provides 
further motivation to consider the supersymmetric extension of the model (SSM+N). Then, 
leptogenesis is induced in both singlet neutrino and singlet sneutrino decays. The results 
are modified by factors of order one, but the basic mechanism and the order of magnitude 
of the asymmetry remain very much the same as in the non-supersymmetric version. 

Supersymmetry must, however, be broken. In addition to the soft supersymmetry break- 
ing terms of the SSM, there are now terms that involve the singlet sneutrinos N, in particular, 
bilinear (B) and trilinear (A) scalar couplings. These terms provide yet another source of 
lepton number violation and of CP violation. One may ask whether these terms can play 



a significant role in leptogenesis. One finds that for a certain range of parameters, the soft 

eptogenesis 



breaking terms play a significant role, and may even be dominant in 
scenario has been termed 'soft leptogenesis.' (For related work, see 

n 

In p] we investigated soft leptogenesis related to CP violation in mixing (a leptonic analog 
of Re(e)x^Tr^)- In this work, we present all the contributions to the lepton asymmetry 
that arise in this scenario. Tne OO—O. io B do— over the other 

contributions in a large part of the parameter space. Yet, if the scale M is relatively low, 
other contributions, related to CP violation in the interference of decays with and without 
mixing (a leptonic analog of Sb^k), and to CP violation in decay (a leptonic analog of 
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Re(e')-K'- + T7r)j P^ a y a significant role. 

The plan of this paper is as follows. In section |H] we derive exact expressions for the 
singlet sneutrino decay rates into final (s)leptons in terms of mixing and decay amplitudes. 
In section Hm we present our model, that is the supersymmetric standard model extended to 
include singlet neutrinos (SSM+N) and express the mixing and decay amplitudes in terms 
of the model parameters. Our main results are obtained in sections IIVI and In section 
IIVI we evaluate the lepton asymmetry in terms of the model parameters and, in particular, 
assuming hierarchy between the supersymmetry breaking scale and the mass scale of the 
singlet sneutrinos, find the potentially leading contributions. In section we estimate the 
size of the various contributions and find the regions in the SSM+N parameter space where 
these contributions can account for the observed baryon asymmetry. We summarize our 
results and draw further conclusions in section IVI1 Additional points are made in two 
appendices. In appendix [X] we explicitly prove that the consideration of three body final 
states does not change the picture. In appendix |B] we discuss the possibility of naturally 
obtaining a small B term for the singlet sneutrinos. 



II. MIXING AND DECAY 



We would like to calculate the CP-violating lepton asymmetry: 

£ ^ r(Z) + r(L)-r(Zt)-r(g 
1 r(L) + r(L) + r(Lt) + r(L)' 

where T(X) is the time-integrated decay rate into a final state with a leptonic content X. 
Here L(L) is the (anti)lepton doublet and L(L') is the (anti)slepton doublet. 
A crucial role in our results is played by the N — A 7 "' mixing amplitude, 

(N\H\Ni) = M 12 - ~r 12 , (2) 

which induces mass and width differences, 

= AM = M H -M L _ AT _ Th-Tl , . 

x- T - r , y - 2r - 2r , [6) 

(r is the average width) between the two mass eigenstates, the heavy \Nh) an d the light 

Wl), 

\N L)H )=p\N)±q\m). (4) 



The ratio q/p depends on the mixing amplitude ratio: 
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q 



2M* 2 - iT\ 2 



2M l2 -iT l2 ' (5) 

For each final state X, we define a pair of amplitudes and a quantity Xx involving the 
amplitude ratio and the mixing amplitudes: 

qA x 



A x = (X\H\N), A x = (X\H\Ni), X 



x 



pA 



x 



(6) 



Defining \N(t)) and \N^(t)) to be the states that evolve from purely \N) and \N^), 
respectively, at time t = 0, we obtain the following time-dependent decay rates into a final 
state X: 



2-rt 



r(N(t)^X) = Af x \A x \ 2 e 



+ TZeXx sinh 
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cosh 
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XmXx sin(AM t) 



T(N\t)-^X) = M X \A 
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e~ rt 
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1 + IA 
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cosh 



AH 1 - |A 



x I 



+ IZeXx sinh 



An 



2 2 
ImXx sin(AM t) 



cos(AM t) 



(7) 



where Mx is a phase space factor. Summing over the initial states, N and W , we obtain 
the following four time-integrated decay rates (in arbitrary units): 



T(L) = KIA, 



(l + |p/g| 2 )(l + |A z | 2 ) , (l-b/g| 2 )(l-|Ar| 2 ) 



2(1 -y 2 ) 



2(1 + x 2 ) 



y(l + \p/q\ 2 )neX z x(l - \p/q\ 2 )lmX } 



1-y 2 



1 + x 2 



r(Lt) = M S \A 



Ltl 



(l + |g/p| 2 )(l + |Ar t r 2 ) , (l-|^| 2 )(l-|Ar t |- 2 ) 



2(1 -y 2 ) 



+ 



2(1 + x 2 ) 



+ 



2/(1 + k/p| 2 )^e^ z(l - |g/p| 2 )Xm^- 

it z/t 



1 + x 2 



r(L) = a/})^ 



(l + b/g| 2 )(l + |A r | 2 ) (l-| p / g p)(l-|A r | 2 ) 
2(1 -y 2 ) 2(1 + x 2 ) 

y(l + \p/q\ 2 )KeX T x(l - \p/g\ 2 )TmX T '] 
1 — y 2 1 + x 2 

(l + \q/p\ 2 )(l + \X L \- 2 ) , (l-\ q /p\ 2 )(l-\X L \- 2 ) 



r(L) = a^i^ 12 



+ 



2(1 -y 2 ) 2(1 + x 2 ) 

+ |g/p| 2 )fte^ z(l - \q/p\ 2 )l mx -' 



1 — y 2 1+x 2 
Using these four decay rates, we can obtain an exact expression for defined in eq. (JTJ) 



(8) 
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III. THE SSM+N 



Since we are interested in the effects of the soft super symmetry breaking couplings, we 
work in a simplified single generation model. The relevant superpotential terms are 

W = Ye Q/3 L a NHp + ^MNN, (9) 

where L is the supermultiplet containing the left-handed lepton doublet fields, N is the 
superfield whose left-handed fermion is the SU(2) x {7(l)-singlet z7 L , and H is the Higgs 
doublet (usually denoted by H2). The relevant soft supersymmetry breaking terms in the 
Lagrangian are the following: 

£ssb = - (m 2 X a 2 X a 2 + Ae a/3 L a NH(i + BNN + h.c.) . (10) 



Here A2 (a = 1,2,3) are the SU(2)l gauginos, N,L,H are scalar fields (and N,L,h are 
their fermionic superpartners). The U(l)y gaugino, Ai, would give effects that are similar 
to those of A2 and can be included in a straightforward way. 

The Lagrangian derived from eqs. (jHJ) and (jlOj) has two independent physical CP violating 
phases: 

<f) N = arg(AMB*Y*), 

<j) W = arg(m 2 M5*). (11) 

These phases give the CP violation that is necessary to dynamically generate a lepton 
asymmetry. If we set the lepton number of and N to —1, so that Y and A are lepton 
number conserving, the two couplings M and B violate lepton number by two units. Thus 
processes that involve Y or A, and M or B, would give the lepton number violation that is 
necessary for leptogenesis. 

There are several dimensionful parameters in Q and (jlOj) . Of these M is supersymmetry 
conserving and all other are supersymmetry breaking. We assume the following hierarchies: 

fsS ^«l, (12) 

where msusY is the supersymmetry breaking scale in the SSM+N (we take msusY ~ 1 TeV), 
and, unless otherwise stated, 

\m 2 \ ~ \A/Y\ ~ \B/M\ ~ m SUSY . (13) 




FIG. 1: Two-body decay diagrams of a singlet sneutrino. 



We also assume that \Y\ <C 1, as is required by the condition of out-of-equilibrium decay 
[see eq. (ETTj)]. 

We can evaluate the various parameters of eq. (jHJ in terms of the Lagrangian parameters 
of eqs. (JHJ) and (fTUJ) . The singlet sneutrino decay width is given, for \MY\ ^> \A\, by 

_ IMF 2 ! 



4tt 



For the mixing parameters, we obtain 
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As concerns the decay amplitudes, CPT guarantees the following relation: 



(14) 



(15) 



\A^+\A l ^ + \A T \" + \A l 



\A 1 \ 2 + \A~ L ^ + \A T \ 2 + W 2 



(16) 



We consider only two body final states, since three body (or higher) states give only small 
corrections, as shown in Appendix EI In Fig. ^ we show the relevant diagrams (including 
the dominant one loop corrections) for the four two body final states: (1) LH, (2) U H\ (3) 



Lh, and (4) Lh. The amplitudes are given by (we use an = sign when the difference between 
the absolute values of two CP conjugate amplitudes is negligible and an ps sign when it is 
not) 

\MY\, 
\MY\, 





= I^Zt 


\A T \ 


~ \A L \ 
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\A T \ 
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- \M 
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|m 2 y|^/A 2 + / 2 2 = <!?(a 2 |m 2 y|) 



m 2 A 



M 

m 2 Y\f 1 sm((j) W 



fi sin(0vy - 4) N ) = 0(a 2 \m 2 A/M\), 
b N ) = O (a 2 \m 2 Y\) , 



(17) 



where a 2 = g 2 /(4ir) is the weak coupling constant and where we define 
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h 



In 
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mo 



ml 



+ Li 2 



(18) 



/ l + 4m|/M 2 / \\ + yfl + 4ml/M 2 / 

These expressions assume, for simplicity, that m 2 > mj-,m^, thus neglecting corrections 
proportional to m^ and m#. The function Li 2 (z) = MlzHfff j s th e dilogarithm function. 
For the relevant strong and weak phases, we obtain 

1 



(f) f 



2 ar S ( A Z A Zt 
- arg (A r A L ) 



sin <5 S = sin ■ 



arg(A z A ?t 1 ) 3« 2 



m 2 Y 



A 



sin 5/ = sin 



arg(A r A L 



,fl 



/1 = e>K 

= 0(1). 



(19) 



7i 2 + fi 

Note that there are several relations between decay amplitudes to final scalars and to final 
fermions. These relations have to be taken into account when evaluating the asymmetry. 
First, in the supersymmetric limit we have \Ar\ = \A-^\. Second, we have \Aj j /Aj\ { sm5 s = 
\Aj i j A-j\ sin5j. Similar relations hold for the CP conjugate amplitudes. 
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IV. THE LEADING CONTRIBUTIONS TO e e 



Many terms that contribute to the lepton asymmetry are small and can be neglected. 
The small parameters that play a role are the ratio eg, the weak coupling constant ct2, and 
the Yukawa coupling Y. The dependence on the Yukawa coupling enters either via the 
combination |A|/|MY| which, as can be seen from eq. (|13|). is taken to be of order es, or 
via the x parameter evaluated in eq. ([15)1 . The x parameter can be small or large but, since 
x ~ 87res/|^| 2 , we take x ^> es- (Some of the contributions that we consider are significant 
only for B <C e$M 2 and, consequently, x <C 87reg/|Y| 2 . In these cases, however, x ~ 1 is 
required, so that x ^> e s is still valid.) We keep the x dependence explicit. 

We identify several interesting contributions to eg. We write down only the potentially 
leading contributions and neglect terms that are suppressed by higher powers of es and/or 
«2- We classify the contributions according to the source of CP violation: 

(i) CP violation in mixing: Here, CP violation comes from \q/p\ 7^ 1 (as in Re(e) in 
K — > ixiv). We identify two potentially significant contributions. The first is given by 

2 „ 2 \ /WV.^N 

(20) 
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xA sf e s 
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(21) 



K{\A~ L 


2 + 






2 + 


A L 


2 ) 


M S {\A~ L 


2 + 




2 )+m^ T 


2 + 


Al 


2 ) 



At zero temperature, A s f = 0(e|), but for temperature at the time of decay that is compa- 
rable to the singlet sneutrino mass, Td ~ M, we have A s / w (Af s — N/)/{N s + A//) = 0(1). 
The second contribution is given by (neglecting now corrections of order A s f) 



m 
5 2 



r 



4 1 



O 



x z 
_e|a2_ 

1 + X 1 



A- 



A-, 



A 



it 



,4 



Lt 



sin 5* cos d>* — 



A T 



At 



A, 



A, 



sin 5f cos 



(22) 



(ii) CP violation in interference of decays with and without mixing: Here, CP violation 
comes from arg(AxA^) 7^ (as in S,i, k s in B — >• J /ipK and similar to the mixing contribution 
to standard leptogenesis (see e.g. |li|), though mixing in the latter case is between different 
generations rather than between CP conjugate states). We identify the following potentially 
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significant contribution: 
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(iii) CP violation in decay: Here, CP violation comes from \Ax\ ^ \A-^\ (as in Re(e') 
in K — > tttt and as in the vertex contribution to standard leptogenesis). We identify the 
following potentially significant contribution: 
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Lt 
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Lt 
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(24) 



(iv) We also find a contribution that involves all three types of CP violation and is not 
necessarily sub-dominant: 
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4(1 +X 2 
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Lt 
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(25) 



We note that, apart from e™, all the contributions involve loop diagrams with gaugino 
exchange. The gaugino is playing a double role here. First, its mass provides a new physical 
CP violating phase. Second, the loop diagrams provide a strong phase. Consequently, direct 
CP violation becomes a possible source of the lepton asymmetry. Gaugino interactions do 
not violate lepton number, but they allow the lepton number violating time evolution of the 
heavy sneutrinos to contribute to ei in new ways. Without gaugino interactions, indirect 



CP violation is the only significant source of soft leptogenesis la, LZJ • Direct CP violation 
can still be induced, but it involves higher powers of the Yukawa couplings and is therefore 
negligibly small. 



V. THE SIZE OF e e 

In the previous section, we distinguished five potentially important contributions to Eg_. 
These five contributions can be separated into three different classes: 

e t = e? + (e< + e d ) + {e™ + e mdi ), 

.d /r\ ( 2 



O (4aa , 



e^e^ = o(jfe)- (26) 
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The generated baryon to entropy ratio is given by 

n B /s ~ -k10~ 3 q, (27) 

where k < 1 is a dilution factor which takes into account the possible inefficiency in the 
production of the heavy sneutrinos or erasure of the generated asymmetry by lepton num- 
ber violating scattering processes. Since observations determine %/s ~ 10~ 10 , any of the 
contributions in ()26j) would be significant only if it yields \si\ > 10~ 6 . We now specify 
the conditions on the parameters whereby each of the three classes of contributions can be 
responsible for a successful leptogenesis. Since all the effects that we consider are related 
to supersymmetry breaking and therefore suppressed by powers of es, soft leptogenesis can 
give significant effects only for es > 10~ 6 , that is, 

M < 10 6 m SUSY ~ 10 9 GeV. (28) 

In order that the singlet neutrino and sneutrino decay out of equilibrium, we should have 
a decay rate, T = M\Y\ 2 /47T, that is not much faster than the expansion rate of the Universe, 
H = 1.66g¥ 2 T 2 /rrip\ (g* counts the effective number of spin degrees of freedom in thermal 
equilibrium; g* = 228.75 in the SSM), at the time when the temperature is of order M: 

M/\Y\ 2 > 3 x 10 16 GeV. (29) 

On the other hand, the sneutrino decay should occur before the electroweak phase transition, 
when sphalerons are still active, T > H(T ~ 100 GeV): 

M\Y\ 2 > 3 x 10~ 13 GeV (30) 

Combining eqs. (|2~%j) . (f2~9"j) and (JBTij) . we learn that soft leptogenesis can give significant 



effects only for 



„ ,VPGeV\ ip - . ( M ^ 



10 *w*^Kwg&) ■ (31) 



With such a small Yukawa coupling, the decay width is rather narrow, 

/ M \ 2 

(i) The contribution from e" 1 is of order (x/(l + x 2 ))A s f€s. For temperatures well below 
the mass M, the finite temperature contribution to A s f is given by the following approxi- 
mation (n sJ = (e M /( 2T ) T l)- 1 ): 

A s/ - { ' + ns \ 2 ~ { '~ nf l - 2e~^\ (33) 
si (1 + n s ) 2 + (1 - n f ) 2 V ; 
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where Td is the temperature at the time of decay. To obtain \sg\ > 10 6 we must have 

T 1 

— > . (34) 

M ~ 21n(2e s /10- 6 ) 1 ; 

By using V = H(Td), this can be translated into an upper bound on M/\Y\ 2 : 

M/\Y\ 2 < 4 x 10 16 Ge1/[21n(2e s /10- 6 )] 2 . (35) 

The lower (|35J) and upper (}2T)|) bounds define, for given M, a range for \Y\ and a range for T. 
Finally, we must have x/(l + x 2 ) > 10~ 6 /(A s yes). Taking into account that x = 2\B\/(MT), 
for a given value of M we obtain an allowed range for B. Since the naive estimate is 
\B\ ~ MmgusY, h is useful to write the allowed range for \B\ in units of Mm^vsY- We do 
so in Fig. 121 We conclude that e™ can account for the observed baryon asymmetry under 
the following conditions: 

1. The mass of the lightest sneutrino is light enough, M < 10 9 GeV. 

2. The Yukawa couplings are small enough, Y < 10~ 4 . The lighter is M, the smaller the 
Yukawa coupling must be. 

3. The B parameter is well below its naive value, \B\/(Mms\jsY) 10~ 3 . The lighter is 
M, the more suppressed the B coupling must be. 

We note that the inclusion of three body decays Q] does not change the basic picture and, 
in particular, does not modify the estimate of A s f. We prove this statement in Appendix 

El 

(ii) The contribution from e d and e % is of order a 2 e| . Since «2 ~ 10 -2 , we must have 

M<10 2 m SU SY. (36) 

Eq. (|3*T]) then requires 

Y < 10" 6 . (37) 

The region where this class of contributions can account for the observed baryon asymmetry 
is to the left of the dash-dotted line in Fig. El Note that \B\ is not constrained in this 
scenario. In particular, it can take its naive value, \B\ ~ MmsusY; in which case x = 
2\B\/(MT) > 10 11 , so that the sneutrino oscillation rate is much faster than its decay rate. 
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FIG. 2: Regions in the Y — B plane where T < H(T = M) and e t > 10~ 6 . We take m SU SY = 10 3 
GeV. The approximation made in our calculations (x > es) does not hold below the dotted line, (i) 
For E£ ~ e™, the allowed regions are within the solid curves, for M = 2 x 10 8 (right), 10 (middle) 
and 10 5 (left) GeV. (ii) For e t ~ e d + e\ the allowed region is to the left of the dash-dotted line, 
with M = 10 5 GeV. (iii) For e e ~ ef + e mdi , the allowed region is below and to the the left of the 
dashed curve, for M = 10 5 GeV. 

(iii) The contribution from and £ md * is of order a2e|/(l +x 2 ). Consequently, the bound 
(jnnj) on M and the bound (|H7jl on F apply. In addition, we must have x ^> 1, which implies 



This third class of contributions is never much larger than the second class. It may, however, 
be comparable if B is small enough. The region where this class of contributions is significant 
is to the left and below the dashed line. 

We note that, since our calculations are performed with the assumption that x ^> es, 
they should not be trusted for B/(Mms\jSY) < Y 2 /(8tt), that is below the dotted line in 



B MY 2 



n 



(38) 



< 



— < 10 



Mmsusv ~ ^susY 8tt 



Fig. El 



VI. CONCLUSIONS 



Our main conclusions regarding the range of parameters where soft leptogenesis may be 
successful are the following: 
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1. Soft leptogenesis can be neglected for M S> 10 9 GeV. 

2. Soft leptogenesis can work for M 3> 10 5 GeV only if the Yukawa couplings have small 
values in a rather narrow range and if the B parameter is very small compared to 
its naive scale (Mmsusy)- (We comment on the possibility of naturally achieving 
B -C MmgusY in Appendix |Bj) 

3. For M < 10 5 GeV there are several contributions from soft leptogenesis that could 
account for the observed baryon asymmetry. All the supersymmetry soft breaking 
terms can assume their natural values. 

The main novel point of this paper is the realization that soft supersymmetry breaking 
terms give contributions to the lepton asymmetry that are related to CP violation in decays 
[e d of eq. ()24|)] and in the interference of decays with and without mixing [e % of eq. ([23)1] . 
In contrast to CP violation in mixing [e m of eq. ()2()j1 ]. the oscillation rate needs not be 
comparable to the decay rate in order to have a significant effect. This is the reason that 
the B term can assume natural values. The new contributions to eg are second order in 
supersymmetry breaking terms and further suppressed by a loop factor (~ a 2 ^2^4/ (M 2 Y)) 
and are, therefore, significant only if M is not much higher than 10 2 msusY- 

The contribution to the lepton asymmetry related to CP violation in mixing (e™ of eq. 
(I2()|) ). which was originally discussed in refs. , requires thermal effects in order to be 

significant. In contrast, the new contributions discussed here (such as e % of eq. (f2*3"|) and e d 
of eq. (j2U)) do not require thermal effects and, consequently, allow a non-thermal scenario of 
leptogenesis to work. Such a scenario would arise if, for example, sneutrinos were produced 
by inflaton decays (or if the sneutrino itself were the inflaton) , and the temperature of the 
thermal bath at the epoch of decay is well below M (though above the electroweak scale so 
that sphalerons are still active). 

Soft leptogenesis opens up the interesting possibility that the scale of the lightest singlet 
(s)neutrino mass (M) is not far above the electroweak scale. In contrast, standard leptoge- 
nesis cannot yield, in general, a large enough asymmetry for low M. The difference between 
the two scenarios lies in the different role of the Yukawa couplings. In both standard and 
soft leptogenesis, the condition for out of equilibrium decay associates a low scale M with 
tiny Yukawa couplings Y. In standard leptogenesis, the Yukawa couplings are the source 
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FIG. 3: Three-body decay diagrams of a singlet sneutrino 

of CP violation; therefore, small Y yield a small e£. In soft leptogenesis, CP violation is 
induced by soft super symmetry breaking terms and is not suppressed by small Y. 

APPENDIX A: THREE BODY DECAYS 



The H u field has Yukawa couplings to neutrinos and to up quarks. The superpotential 
terms, W = YNLH + Y u QuH, give a quartic scalar interaction term in the Lagrangian, 



£ 4 = YY*NLQ^v) + h.c. 



(Al) 



where Q is the scalar quark doublet and u is the up-singlet. This coupling allows the three 
body decay mode, N — > L'uQ. Since there is no similar quartic coupling of N to two 

(3) 

fermions and one scalar, one may think that for the three body decays, the vanishing of A^y 
(defined in (|A3|) ) in the supersymmetric limit is avoided, and a sizeable lepton asymmetry 
is induced even at zero temperature. This is, however, not the case, as we now explain. 
We are considering contributions to the CP asymmetry of the form 



r (3) 



X 



4(1 



a(3) 



where 



A {3) = 



S? =1 (-i) L ' + W<|4 3) l 2 
Af s \A L \ 2 +Af f \A T \* 



(A2) 



(A3) 
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Here A\ is the amplitude of a relevant final three body state with lepton number Li = ±1. 
At zero temperature, all the three body phase space factors Mi are equal in the supersym- 
metric limit and, consequently, A^j oc J2t=i(~ l) ii+1 |4 3 ' ) | 2 - The five tree level diagrams, 
leading to four different final states, are shown in Fig. El A straightforward calculation 
gives: 



4 



(3) 



A (3) + 4 (3) 



2\Y\ 2 \Y,. 



m 



23 



23 



2^2 : 



(3) 



.4 



4 3) 



A 



(3) 



= 2\Y\ 2 


Y u \ 


= 2\Y\ 2 


Y u \ 


= 2\Y\ 2 


Y u \ 



12" '23 



(m 2 23 - y 2 ) 2 " 



2 2 
m l3 m 23 

(mh-^r 
M 2 m\ 3 



[m 



23 



2\2 ' 



(A4) 

(fcj + kj) 2 , with 



where /x is the supersymmetric mass of the H u supermultiplet, and mfj 
ki,k 2 ,k 3 the momenta of, respectively, the final (s)lepton, the singlet up (s)quark and the 
doublet (s) quark. Then, 



E(-i) Li+1 l4 

i=l 



(3) 1 2 



l4 3) l 2 + 

2\Y\ 2 \Y,. 



i4 3) r+i4 3) r 



L4 



(3) 1 2 



2 m 2 23 (M 2 - mj 2 - m 2 3 - m 2 23 ) 



(m 2 23 - ft 2 ) 2 



0. 



(A5) 



The last equation, that is the vanishing of the 



4 ( 



\Li+l\ 4(3)|2 



A\ I , holds in the super- 



symmetric limit, when the three final particles are massless. The result is that, in the 

(3) (3) 
supersymmetric limit, e\ — 0. The vanishing of A^j is lifted by finite temperature effects, 

similarly to the case of A s /, but then the contribution of the three body states is small 

compared to the dominant two body ones. 

If we assign lepton number L = to the 7V-supermultiplet, then the quantities A s f 

defined in eq. 023J and A J } defined in eq. f|A3jl are the asymmetries between AL = +1 

and AL = — 1 decay rates. Then, the vanishing in the supersymmetric limit of A s f and 

(3) 

A s j , demonstrated explicitly in our work, becomes understandable on general grounds and 
generalizes to n-body states for any n. In a single generation framework and in the absence 
of supersymmetry breaking, singlet neutrino decay rates to leptons and antileptons must be 
equal. Then, by supersymmetry, this should hold also for singlet sneutrinos. 
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APPENDIX B: ON THE NATURALNESS OF B = 



We consider the following superpotential terms: 

W = MNN + YNLH, (Bl) 

and SUSY breaking terms, 

C = BNN + ANLH. (B2) 

In the absence of these terms, there are four additional flavor conserving global U(l) sym- 
metries: t/(l)jy x f/(l) L x U{1)h x t/(l)_R, with the following charge assignments: 

#(1,0,0,0), #(1,0,0,1), 
£(0,1,0,0), £(0,1,0,1), 

/i(0, 0,1,0), H{0, 0,1,1). (B3) 

Selection rules for the symmetries may be used if M, Y, A and B are treated as spurions 
with charges assigned to compensate those of the fields: 

M(-2, 0,0,0), F(-l, -1,-1,-1), ^(-1,-1,-1,-3), £(-2,0,0,-2). (B4) 

To understand the consequences, it is simpler to examine the charges of the spurions under 

U(1) N - L X U(1) 2 R-3(L+H) X U(1) 2R -(L+H) X U(l) L+H \ 

M(-2, 0,0,0), F(0,+4,0,0), A(0, 0,-4,0), B(-2, -4, -4, 0). (B5) 

We learn that setting 5 = does not add a symmetry to the Lagrangian. Consequently, B 
is additively renormalized. However, setting B and any other of the three couplings, M, Y 
or A, to zero is natural. 

We can therefore think of a three generation framework where, for example, Y = 



because of a supersymmetric Froggatt-Nielsen symmetry 



15, UM. Then = is natural. 



When the FN symmetry is spontaneously broken, B can be naturally suppressed: 

B oc AMY ] < AM. (B6) 

Of course, a Froggatt-Nielsen symmetry can also induce A ~ m SU sYY ?^susy, leading to 
further suppression of S compared to MmsusY- Both the additive renormalization, and the 
suppression factors in (|B6|) are manifest in the RGE [l7|: 

16tt 2 — 5 = MY* A, (B7) 
16 



If, however, B is radiatively generated, as in (|B7J) . the phase <pN vanishes at this order. At 
two loops, there will be a contribution to B that depends on rri2, but then <pN ~ 
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